Fermi liquid in the Hubbard Model with an electron reservoir: 
Normal state of cuprate superconductors 
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It is proved that the ground state under the supreme single-site approximation (S'^A), the dy- 
namical mean-field theory (DMFT), or the dynamical coherent potential approximation (DCPA) is 
the normal Fermi liquid in the presence of an infinitesimally weak hybridization with an electron 
reservoir, except for the just half filling of electrons and the infinite on-site repulsion. In the strong- 
coupling regime, in particular, the Fermi liquid is stabilized under S'^A, DMFT, or DCPA by the 
Kondo effect, which stabilizes a local singlet on each unit cell, and is further stabilized beyond it by 
the Fock-type term of the superexchange interaction or a resonating valence bond (RVB) mecha- 
nism, which stabilizes a local singlet on each pair of nearest neighbors. The Fermi liquid is a relevant 
normal state to study possible lower-temperature phases or the true ground state. It is proposed 
that the Fermi liquid stabilized by the Kondo effect and the RVB mechanism is the normal state of 
cuprate high-temperature superconductors. 



PACS numbers: 71.10.- 



71.10.Ay, 71.27.-fa, 74.20.- 



I. INTRODUCTION 

High temperature (high-Tc) superconductivity in 
cuprates is an interesting and important issue in solid- 
state physics J^i^i^i^i^ Parent cuprates are insulators, 
which show antiferromagnetism at low temperatures. 
When holes or electrons are doped, they become metals, 
which show exotic properties i^'^'^i^ High-Tc supercon- 
ductivity occurs in exotic metals in the vicinity of the 
Mott metal-insulator (M-I) transition, which is also an 
interesting and important issue i^i^iSii^ In order to resolve 
the issue on the mechanism of high-Tc superconductivity, 
it should be clarified whether the normal state above Tc 
is an exotic Fermi liquid (FL), such as the resonating va- 
lence bond (RVB) statC)^ or the conventional or normal 
FL. First of all, the nature of electron correlations in the 
vicinity of the Mott transition should be clarified. 

The Hubbard model is an effective Hamiltonian for 
the Mott transition. According to Hubbard's theory,— li^ 
when [/ > IF, with U the on-site repulsion and W the 
band- width, the band splits into two subbands, i.e., the 
Hubbard gap opens between the upper Hubbard band 
(UHB) and the lower Hubbard band (LHB). When n = 1, 
with n the electron density per unit cell, the ground state 
is a prototype of the Mott insulator, which seems to be 
an abnormal insulator characterized by the ground-state 
entropy diverging in the thermodynamic limit j^ the in- 
sulator for n = 1 and U /W — +oo is the typical Mott 
insulator, whose entropy is /cb In 2 per unit cell. When 
[/ < IF or n 7^ 1, the ground state is expected to be a 
metal since the density of states (DOS) at the chemical 
potential is nonzero; the Fermi surface (FS) cannot be 
defined within Hubbard's theory. 

According to Gutzwiller's theoryji^ii^ii^ together with 
the FL theory, ^'*''-^- the FS is defined in the quasi-particlc 
band, which is called the Gutzwiller band in this pa- 
per. According to Brinkman and Rice's theory^i^ how- 
ever, when n = 1 the Gutzwiller band vanishes and an 



M-I transition occurs aX U — C/br, with C/br — W . The 
specific heat coefficient is diverging as C/ ^ C/br — 0, 
which implies that the ground state is the Mott insulator 
or the abnormal insulator when n = 1 and U > C/br- 
The ground state is the FL when n y^ 1 or U < C/br- 

When U is large but is still finite, the perturbative 
process that gives the superexchange interaction never 
vanishes^ so that matrix elements must be nonzero 
among the degenerate ground states if their configura- 
tions are close to each other. It should be critically ex- 
amined whether the abnormal insulator is really stable 
for finite U or the third law of thermodynamic is really 
broken for finite U, but within the restricted Hilbert sub- 
space where no order parameter exists; it is obvious that 
the abnormal insulator for finite U is unstable against an 
antifcrromagnetic state in the whole Hilbert space. 

One may speculate that DOS has a three-peak struc- 
ture, with the Gutzwiller band between UHB and LHB, 
in a metallic phase of ti ~ 1 and U > W. Hubbard's and 
Gutzwiller's theories are under the single-site approxi- 
mation (SSA). According to another SSA theoryf2i the 
Gutzwiller band appears at the top of LHB for n < 1, 
which implies that it appears at the bottom of UHB for 
n > 1. The SSA that considers all the single-site terms 
is rigorous for d -^ +oo within the restricted Hilbert 
subspacer2i22ii2ii2^i^ with d being the spatial dimension- 
ality. The SSA is called the supreme single-site approxi- 
mation (S^A) in this paper. The S^A is reduced to solv- 
ing the Anderson modelf^li^i^Sii^S which is an effective 
Hamiltonian for the Kondo effect. The three-peak struc- 
ture corresponds to that in the Anderson model, with the 
Kondo peak between two subpeaks. The Kondo effect 
has relevance to electron correlations in the vicinity of 
the Mott transition. The S'^A is also formulated as the 
dynamical mean-field theory^i^ai^i^^ (DMFT) and the 
dynamical coherent potential approximation^ (DCPA). 

A Kondo-lattice theory21i^i2£ and a cluster DMFT 
(CDMFT>2ii^i26i37 j^j.g proposed beyond S^A or DMFT. 



In the Kondo-lattice theory, it is assumed that the ground 
state under S^A or DMFT is the FL, i.e., the Gutzwiller 
band and the FS exist or they survive even if the Hubbard 
gap opensi^i^^i^^ Numerical resuhs of S'^ A or DMFT and 
those of CDMFT give an indication that, when n ~ 1 and 
U/W > 1, the Gutzwiller band vanishes and the ground 
state is an insulator ^ I'^^i^^i'^'^i^^i'^^i^'^' i.e., the Mott insu- 
lator or a spin liquid. It is a crucial issue which is the 
ground state under S^A, DMFT, or DCPA, or within the 
restricted Hilbert subspace, the normal FL, an exotic FL, 
the Mott insulator, or a spin liquid. This issue is related 
to that on the normal state of cuprate superconductors. 

The t-J model is an effective Hamiltonians for cuprate 
superconductors, which are anisotropic quasi-two dimen- 
sional oxides composed of CUO2 planes. It is derived 
from the Hubbard model^ and the d-p model^^ which 
considers d orbits on Cu ions and p orbits on O ions on 
Cu02 planes. When n = 1, the t-J model is reduced 
to the Heisenberg model. Anderson proposes the RVB 
theory of high-Tc superconductivity in the t-J modelj^ 
A parent state is a spin liquid or the RVB state*° in the 
Heisenberg model, rather than the Mott insulator. When 
holes or electrons are doped, it becomes a metallic RVB 
state, which is the normal state in the RVB theoryii^ 
The RVB state, which is an insulator or a metal, is sta- 
bilized by the formation of a resonating valence bond 
due to the superexchange interaction. In the mean-field 
RVB theory)^ the metallic RVB state is stabilized by 
the Fock-type exchange interaction. According to the 
Kondo-lattice theory for the t-J model^^ the FL is sta- 
bihzed under S^A, DMFT, DCPA by the Kondo effect 
and is further stabilized beyond it by the Fock-type ex- 
change interaction. The RVB state and the FL are sim- 
ilar to each other. Similarities and differences between 
the RVB state and the FL should be clarified. 

This paper is organized as follows: In Sec. |lll it is 
proved that the ground state of the Hubbard model un- 
der S^A, DMFT, or DCPA is the FL when an infinites- 
imally weak hybridization with an electron reservoir ex- 
ists, except for ri = 1 and U/W — -foo. In Sec. IIIIl 
it is shown that the FL, which is further stabilized by 
the Fock-type exchange interaction beyond S'^A, DMFT, 
or DCPA, is eventually unstable against at least a mag- 
netic or superconducting state in two dimensions and 
higher. In Sec. |TVl several issues are discussed: impossi- 
bility of the Mott insulator for finite U, the normal state 
of cuprate superconductors, and so on. Conclusion is 
given in Sec. |Vl An inequality, which is used in the proof 
in Sec. [TTl is proved in Appendix |^ In Appendix [B] 
the FL theory for a spin liquid in the Heisenberg model 
is developed and it is proposed that the insulating and 
metallic RVB states are the spin liquid and the FL, re- 
spectively. In Appendix [Cl on the basis of an analysis 
that single-particle excitations are different between the 
presence and absence of an electron reservoir, it is shown 
that the well-known physical picture for the Mott transi- 
tion, which is one in the absence of an electron reservoir, 
is never relevant to explain the Mott transition. 



II. FERMI LIQUID UNDER S^A OR DMFT 
A. Perturbation from an electron reservoir 

In this paper, a reservoir is explicitly considered: 

n = H + Hros + \v-^J^r. (2.1) 

The first term is the Hubbard model defined by 

H = ea^Uia + ^ tija\^aj„ + U'^Ui^riii, (2.2) 

with riia = al^Uia, e^ the band center, tij transfer inte- 
grals, and U the on-site repulsion. The dispersion rela- 
tion of electrons is defined by 

E{k) = ea + ^ ^ i^J exp [zk • (R, - R,)] , (2.3) 



«7^j 



with A^c the number of unit cells and R,; the position of 
the ith unit cell. DOS is defined by 



D{e) = ^Y.^l^-^(^)^- 



(2.4) 



The band-width of £:(k) or D{e) is denoted by W. The 
electron density n per unit cell in the Hubbard model is 
defined by n = (N) /N^, with 



AA.^. 



(2.5) 



The chemical potential for C/ = and T — K, which is 
denoted by iJ-a{n), is defined by 



= irE^([^o(")-^(k)]/W^), 



(2.6) 



with 6{x) the Heaviside function, i.e., 6{x) = for x < 
and 9{x) = 1 for a; > 0. It is assumed that D{e) is 
nonzero, continuous, and finite at least at £ = /uo(n). 
The second term stands for the reservoir: 



Kcs-^^6(k)5[^6ka 



(2.7) 



kCT 



The third term is a hybridization term between the Hub- 
bard model and the reservoir: 



AV = A ^ \v. 



4A^ + '"*b]„ai. 



(2.8) 



with the summation over i being over TZ of randomly 
distributed hybridization sites, and 






(2.9) 



It is assumed that A is a nonzero but infinitesimally small 
numerical constant, which is denoted by A = ±0. Unless 



A = 0, the electron number M in the Hubbard model 
is a non-conserved quantity. It is assumed that (( Uj )) — 
((<)) = and ((wjw*)) = %nh|wp, with ((• • •)) standing 
for the ensemble average for TZ and nh the density of 
hybridization sites per unit cell. In the last term, fj, is 
the chemical potential, and 



^ = J2 ('^l"- + ^l^- 



(2.10) 



The Green function for electrons in the Hubbard model 
with U — averaged over the ensemble is given by 



Gi°)(i£„,k) 



1 



ien + ^J.- E(k) + iX^T{ien) ' 



(2.11) 



When A = ±0, the second-order perturbation is accurate 
enough to treat scatterings from the random hybridiza- 
tion, so that 



r(ie„) = inh\v\'^jr-Yl 



1 



N^^ ien + ^^-Eb{k)^ 



(2.12) 



It is assumed that no gap opens in the reservoir or that 
r(e + iO) is continuous at e = and 



Rer(+iO) > 0. 



(2.13) 



For the sake of simplicity, it is assumed throughout 
of this paper that the model (|2.ip is on a lattice in two 
dimensions and higher and tij is only nonzero between 
nearest neighbors (ij) , if nothing is mentioned about di- 
mensionality d or tij; tuj\ is denoted by i, if necessary. 



B. Proof of the FL ground state 

1. Fermi-surface (FS) condition in the Anderson model 

The s-d model is another effective Hamiltonian for 
the Kondo effect. According to Yosida's perturbation 
theory^ and Wilson's renormalization-group theoryjii 
when DOS of the conduction band is nonzero at the 
chemical potential, the ground state is a singlet or the 
normal FL but is exceptionally a doublet for Jsd ~ 0, 
with Jsd the s-d exchange interaction. 

The s-d model is derived from the Anderson model: 

(7 

^ E (^kcLd. + V^dWa) , (2.14) 



kCT 



kcr 



with rida = d\d„^ ed the d electron level, U the on-site re- 
pulsion, Vic the hybridization matrix, and A'a the number 
of unit cells. The hybridization energy is defined by 



^{e) = ^J2\^^\'^[^ + f^-E,{k)], 



Na 



(2.15) 



with fj, the chemical potential. The FS of conduction 
electrons, which is defined by -Ec(k) = fl, exists when 



A(0) > 0, 



(2.16) 



which is called the FS condition in this paper. The results 
for the s-d model imply that, when the FS condition is 
satisfied, the ground state is a singlet or the FL but is 
exceptionally a doublet for n^ = 1 and U/ttA{0) = -|-oo, 
with Ud = '^^ {nda) being the density of d electrons. 

According to Bethe-ansatz solutions for the s-d model 
with a constant DOS of the conduction band and the 
Anderson model with a constant A(e), the ground 
state of either model is the FL except for each ex- 
ceptional casej^i^i^i^ In general, the nature of the 
ground state depends only on relevant low-energy proper- 
ties, such as A(0), and high-energy properties renormal- 
ize only quantitatively the ground state, as is demon- 
strated by renormalization-group theories for the s-d 
model^^i^ The Kondo effect is almost or practically 
solved i^i^i^i^i^ The most fundamental assumption of 
this paper is that when the FS condition is satisfied the 
ground state of the Anderson model is the FL except for 
nd = l and i7/7rA(0) = +oo. 



2. Mapping to the Anderson model 

Consider the Hubbard model within the restricted 
Hilbert subspace where no order parameter exists. The 
Green function averaged over the ensemble is given by 



Gcr{i£n,'^) 



iSn + fJ-- E{k) - T,a{ie„, k) + iX'^r{ien) ' 

(2.17) 
Here, S(j(i£„,k) is the self-energy, which is divided into 
single-site I]o.(Je„) and multi-site A'E^{ien,]i): 



So-(ie„, k) = Scr(ien) + AI]o-(i£„, k) 



(2.18) 



Since any vertex correction due to AV and Tiros is 0(A'*) 
or higher, it can be ignored when A — ±0. The single-site 
S(j(Jen) is given by that of the Anderson model provided 
that the on-site repulsion line and the single-site electron 
lines in Feynman diagrams are the same as each other be- 
tween the Hubbard and Anderson models. If "Ea-ii^n) and 
AI]o-(J£n,k) are obtained, the single-site Green function 
of the Hubbard model is given by 

i?^(«£n) ^^Y. G-(«en, k), (2.19) 



and that of the Anderson model is given by 

1 



G(j(ie„) 



ien + ft- ed- So.(ie„) - - j de'— ^^ 



TT / lEn — e' 

(2.20) 



The condition for the electron Hnes is given by i?o-(«£n) = 
Ga{iSn) or 

R„{e + iO) = Ga{e + tO), (2.21) 

which can never be satisfied unless 

p. - ed = IJL - ta- (2.22a) 

It follows from ((2?2T|l that 

A(e) = Im ta{e + iQ) + R-^{e + iQ) . (2.22b) 

The condition for the on-site repulsion line is given by 

U ^U. (2.22c) 

A problem of calculating the single-site 'E,„{ien) is re- 
duced to that of determining and solving self-consistently 
the Anderson model to satisfy Eq. (|2.22p , which is called 
the mapping condition in this paper; the multi-site 
AEo-(i£n,k) should also be self-consistently calculated 
with the single-site S(j(je„) to satisfy Eq. (|2.22p . Ac- 
cording to Eq. (|2.2ip . DOS and the electron density are 
the same as each other between the two models: 



Pie) 
and 



1 



lniR^{e + iO) 



1 



lmG^{e + iO), (2.23) 



nd 



+ 00 



2 / defie)pis), (2.24) 

J —oc 

l/[exp(e/fcBT) + 1] being the Fermi-Dirac 



with /(e) 
function. 

In S^A, DMFT, or DCPA, only the single-site E^(ie„) 
is considered, so that 



R^{e + iO) = de'- 



D{e') 



(2.25) 
The mapping condition (j2.22bp is iteratively treated to 
determine the Anderson model to be solved. Even if any 
l;^(e -I- iO) is assumed in Eq. (|2.22b|) .^° it follows that 



A(e) >ReA^r(e + iO), 



(2.26) 



as is proved in Appendix |^ According to Eqs. (|2.13[) 
and (|2.26p . the FS condition (|2.16p is satisfied in each it- 
erative process. Therefore, any self-consistent solution is 
of the normal FL. When A = ±0 or A/" is a non-conserved 
quantity, the ground state under S^A, DMFT, or DCPA, 
is the normal FL except for U /W = +oo and n = 1. 



C. Fermi-liquid (FL) relation 

Consider the mapped Anderson model in the presence 
of an infinitesimally small Zeeman energy hz = gfJ-BH 



and an infinitesimally small chemical potential shift Afi. 
The self-energy is expanded at T = K such that 

t,{e + iO) - Eo(0) + (l - 0T,) e + (l - 4>s) ^ahz 
+ (l - <^c) A/i + O (e^) , (2.27) 

with So(0), (j)j, (j)s and cj)c all being real and finite. Since 
DOS of the Anderson model is the same as that of the 
Hubbard model, it follows that 



pio) = ^J2^[^'-^o{Q)-E{■k)]. 



(2.28) 



Physical properties of the Anderson model can be de- 
scribed by the FL relation. ^^ In general, 20^ = (j)s + (j)c, 
so that 1 < (j)s/<P'y < 2; (t)s/4>-y is simply the Wilson 
ratio.i^^^ When n ^ 1 and U jW -^ -|-oo, in particu- 
lar, charge fluctuations are totally suppressed, so that 
(j)c/4>'y — > and 4>s/(t>'y -^ 2. In this paper, the Kondo 
temperature is defined by 



kBTK = l/[4.4>^p{0)]. 
The specific heat coefficient is given by 



2 T.,7 ,„^ TT^fcB 



7 = |7r2fci07P(O) 



6Tk 



(2.29) 



(2.30) 



Physical properties of the Hubbard model can also be 
described by the FL relationiiiii^ According to the FS 
sum rule, 

'^^^ J-E^([/^-^o(0)-^(l^)]/^)- (2.31) 



According to Eqs. ([^ and ((OT|) . 

M-So(0) = /io(n). (2.32) 

According to Eqs. (I23D, fTIE^i and ([O^ . 

p(0) = D [^lo{n)] . (2.33) 

According to Eqs. (|2.22b|) . and (E^TD, 

7rp(0) 



A(0) = 



[ReR„{+iO)f + [np{0)f' 



(2.34) 



It should be noted that neither of p{0) and A(0) depends 
on U when n is kept constant. The specific heat coef- 
ficient is given by Eq. (|2.30p . The Kondo temperature 
Tk, which defined by Eq. (|2.29p . is an energy scale of the 
effective Fermi energy of quasi-particles. 

In the strong coupling regime defined by U > W, DOS 
has the three-peak structure, with the Gutzwiller band 
between UHB and LHB. The Green function ((2T7| is 
approximately given at T = K by 



Gcr(ie„,k) = y- 



ie„-$o(k) 



[incoherent term]. (2.35) 



The first term is tlie coherent term, which describes the 
GutzwiUer band. The dispersion relation and the band- 
width are given by 



eo(k) = [E{k) + So(0) - fi] /^^, (2.36) 



and 



W* = W/(j)^ ~ 4:kBTK. (2.37) 

The incoherent term describes UHB and LHB. 



that J = -it'^/UM According to field theory, the 
superexchange interaction arises from the virtual ex- 
change of a pair excitation of electrons between UHB 
and LHB. 52,53,54 ^j^g^ ^-^^ widths of UHB and LHB are 
considered, \J\ is about a half of it^/U when U ~ W^ 
The second term Jq [iloi , q) arises from the virtual ex- 
change of a pair excitation of quasi-particles. When the 
single-site irreducible three-point vertex function in spin 
channels is denoted by As(ze„7*£n + i^i',i^i), it follows 
that 



III. RVB MECHANISM AND ORDERED 
STATES BEYOND S'^A OR DMFT 

A. Kondo-lattice theory 

The strong-coupling regime is mainly studied in this 
section. The irreducible spin polarization function 
'Ks{iLUi,c^ is also divided into single-site TTs{i^i) and 
multi-site A7rs(ia;/, q): 

TTs{iuiuq) = TTsiitoi) + ATTs{iu!i,q). (3.1) 

The single-site TTsiitoi) is given by that of the Anderson 
model. The spin susceptibilities of the Anderson and 
Hubbard models are given, respectively, by 



and 



Xs{iuJi) 



Xs(JWi,q) 



2Tr sjiuji) 
1 - UiTsiiuJi)' 

2'Ks{iuJi,q) 
1- UTTs{iuJi,q)' 



(3.2) 



(3.3) 



A physical picture for Kondo lattices is that local spin 
fluctuations on different sites interact by an intersite ex- 
change interaction. According to this picture, an intersite 
exchange interaction Ig (ilui , q) is defined by 



Xsiii^m) 



Xs{iuJi) 



1- jlsiiuji,q)xs{i^i) 



(3.4) 



It follows from Eqs. ([STj), ([33), ^^, and dM]) that 

Is{iu;i,q)^2U^A7rAiL0i,q){l + O[l/Uxs{i^i)]}- (3-5) 

When U/W > 1, terms of 0[1/Uxs{i^i)] can be ignored. 
The exchange interaction /s(iw;,q) is composed of 
three terms i^!^ 

Is{iuji,q) ^ Js{q) + jQ{iuJi,q) - 4:A{iuJi,q). (3.6) 

The first term Js(q) is the superexchange interaction: 






(3.7) 



According to the second-order perturbation in t, where 
the widths of UHB and LHB are ignored, it follows 



A,(0,0;0) = (/.,[! -C/^,(0)] 
24>, 



uuo) 



{1 + O[l/C/X.(0)]}, (3.8) 



according to the Ward relation;^ terms of O[l/Uxs{0)] 
are also ignored. When only the coherent part of the 
Green function is considered and Eq. p.Sp is approxi- 
mately used for low-energy dynamical processes, 

jQ{iuJi,q) ^ P{iuji,c[) - —^P{iuji,q), (3.9) 



with 



xm 



{<i>s/4>' 



-E 



/[eo(k)]-/[eo(k + q)] 



(3.10) 



which is derived in the random-phase approximation 
(RPA) for pair excitations of quasi-particles. In Eq. (|3.9p . 
the single-site term is subtracted. The third term 
—4:A{iLUi, q) is the mode-mode coupling term among var- 
ious types of fluctuations. It corresponds to that in 
the self-consistent renormalization (SCR) theory of spin 
fluctuations )^^'^^ which is relevant in the weak-coupling 
regime defined by U/W < 1. 

When Eq. (|3.8p is approximately used, the mutual in- 
teraction mediated by spin fluctuations is given by 

^[U~Xs{0,0;0)f[Xs{itoi,q) - Xsi^LOi)] = ^^^j:{tLOi,q), 

(3.11) 
with 



isii^hqi) 



Is{iuji,q) 
1 - jIs{iuJi,q)Xs{m)' 



(3.12) 



In Eq. p. lip , the single-site term is subtracted and two 
(j)s appear as effective three-point vertex functions. The 
mutual interaction mediated by spin fluctuations is sim- 
ply the exchange interaction /* {iuji , q) . Multi-site or in- 
tersite terms are perturbatively considered in terms of 
Is (iuji , q) or /* [iuji , q) . The perturbative theory is sim- 
ply the Kondo-lattice theory. 



B. Stabilization of the FL by an RVB mechanism 

There are two linear terms in the superexchange in- 
teraction Js{ci), which is the main term of /*(«a;;,q): 
Hartree-type^ and Fock-type'*- terms. The Hartree-type 
term, which gives magnetic Weiss mean fields, vanishes 
if no magnetic order parameter exists; magnetic instabil- 
ities are studied in Sec. IIIICI When only the coherent 
term of the Green function is considered, the Fock-type 
term is given by 



^ / np 



oiE^O^ 



J.(k-p) 



2C7) 



Up)' 



(3.13) 



The factor 3 appears because of three spin channels. The 
Fock-type term AI]o.(k) should be self-consistently cal- 
culated with the single-site I](j(J£n) to satisfy the map- 
ping condition (j2.22bp . The self-consistent S(j(je„) is 
expanded, as it is in Eq. (j2.27p . but with renormalized 
So(0), (/>7, 4>sj and 4>c by the Fock-type term, all of which 
are real and finite.— Then, the dispersion relation of 
quasi-particles is given by 

e(k) = [E{k) + So(0) + AS,(k) - /i] /4>^, (3.14) 

DOS at ^ is given by 

^(0) ^ W E 4^(k) + So(0) + AE,(k) - ^J] , (3.15) 



and the specific heat coefHcient 7 is given by Eq. (|2.30p . 
Consider the two dimensional square lattice {d = 
2). The superexchange interaction is given by Js(q) = 
2 J [cos {qxa) + cos {qytt)], with a the lattice constant. It 
follows that 

AI]o-(k) — -(j)jCjJ[cos{kj;a) -f cos(fcj^a)] , (3.16) 

with 




-eo(k) 



W 



[cos(fca;a) + cos{kya)] . 



Since J is antiferromagnetic, the sign of cj, which is 0(1), 
is such that ASo-(k) enhances the band- width of quasi- 
particles. The Fock-type term depends on d and lattice 
structure. In general, DOS is given by 

p{0)^l/[W + 4>^\cjJ\], (3.18) 

and the band-width of quasi-particles is given by 

W* ~ l/[4>^p{0)] c^ WT^^ + \cjJ\. (3.19) 

The Kondo temperature is given by fceTk — VK*/4. 

Under S^A or DMFT, the FL is stabiHzed by the 
Kondo effect, which stabilizes a local singlet on each unit 
cell. The band- width of quasi-particles is Wjf^^, which 



may be infinitesimally small but is nonzero. Beyond S'^A 
or DMFT, the FL is further stabilized by the Fock-type 
term, which stabilizes a local singlet on each pair of near- 
est neighbors; the band- width is broadened by |cjj|. In 
this paper, the stabilization mechanism is called an RVB 
mechanism, since it is the same as or at least similar 
to that in the RVB theory— If 0^ -^ -|-oo, in particu- 
lar, the FL is totally stabilized by the RVB mechanism 
so that W* ~ \cjJ\. Since DOS is vanishing such that 
PFp(O) —> 0, the ground state is almost a spin liquid, 
which is studied in Appendix iBl but is still the FL, pro- 
vided that U jW is finite. 

It is an issue whether the FL is stable within the re- 
stricted Hilbert subspace where no order parameter ex- 
ists. In addition to the Fock-type term, a few terms for 
the multi-site self-energy, AI](j(e-l-JO, k), are examined 
under an assumption that (/)-y < +00 and no order pa- 
rameter appears. In one dimension, terms proportional 
to £ln|e| appear at T = K,^^ which means that the 
FL is unstable against, al least, an exotic metal or the 
Tomonaga-Luttinger liquid j^^i^'^i^^ In d > 2 dimensions, 
no such term appear, which implies that the FL is stable, 
but within the restricted Hilbert subspace. 



C. Possible ordered states 

Since low-energy excitations are so accumulated in the 
FL that the specific heat is proportional to T at T ^ Tk, 
it is presumably the truth that, in the whole Hilbert 
space, the FL is eventually unstable against an ordered 
state in d > 2 dimensions. When an order parameter 
is specified, it is straightforward to study the response 
function corresponding to it, from which the instability 
condition of the FL against the ordered state can be de- 
rived. In this paper, the phase diagram is out of scope 
but only possible ordered states are examined. 

The magnetic susceptibility Xs(*^ijQ) is given by 
Eq. dnH). Provided that 



[l-i/,(0,Q)x.(0)]^^„j,>0, 



(3.20) 



for any q, the FL is stable against any magnetic state. 
Assume that /s(0, q) is maximal at q = Q. When 
/s(0, Q) or Xs(0) is so large that 



[l-i/s(0,Q)x.(0)]^^^ =0, 



(3.21) 



the FL is unstable below Tn against a magnetic state 
or a spin density wave (SDW) state with Q. Accord- 
ing to the FL relation, ^^ together with Eq. (|2.29p . it 
follows that Xs(0) — I/ZcbTk- In the limit of n ^ 1 
and U jW -^ -l-oo, /cbTk is vanishing so that Xi(0) is 
diverging. The superexchange interaction Js(q) is anti- 
ferromagnetic. When n ~ 1, Jq (0,q) is also antiferro- 
magnetic, as is discussed in Sec. IIVCI When n ~ 1 and 
V jW is large enough, therefore, it is probable that the 
FL is unstable against an antiferromagnetic state. In two 
dimensions, however, Tn = +0 K when the mode-mode 



coupling term, — 4A(0,q), or magnetic critical fluctua- 
tions are self-consistently treated.— 

An average of /^(w + iO, q) over a low-energy region 
|w| < fceTk is expanded as 



Re(/:(u; + zO,q))^^^^/*e^'»-(^-^^ 



(3.22) 



with (• • • )^ standing for the average. Most possible order 
parameters, in addition to a magnetic one, are given by 
the decoupling of 









Q/3 



).(i<T^^)aLa.0a]^a,,, (3.23) 



with <T = {ax,(Ty,az) being the Pauli matrix: (a]^a]^,) 
or (ajr>air) of superconductivity, 'Yliri'^T'^jr) of charge 
bond wave (CBW), and ^^^tt' '^V {o^ir^jT') of spin bond 
wave (SBW), in addition to YIitt' ^V Wir^ir') of niag- 
netism or SDW. 

For the sake of simplicity, the nearest-neighbor com- 
ponent of /*,-, which is denoted by /j*, is only considered. 
When n ~ 1, Js{q) and jQ(0,q) are antiferromagnetic, 
as is discussed above, so that /^ is antiferromagnetic such 
that I^ < 0. When n ~ or n ~ 2, Jq(0, q) can be fer- 
romagnetic, as is discussed in Sec. IIV (Jl When Js(q) is 
weak, I^ can be ferromagnetic such that /j* > 0. Possi- 
ble symmetries or waves of ordered states depend on rf, 
lattice structure, n, the sign of Jj*, and others. 

When II is weak or strong, the FL is unstable against 
an anisotropic superconducting (SC) state, at least, if no 
disorder exists: a singlet one for /^ < and a triplet one 
for Ii > 0. For example, consider the two-dimensional 
square lattice. When n ~ 1, /^ < 0, as is discussed above. 
Then, two singlet waves are possible: anisotropic s wave 
and d'y wave. When any pair breaking by SC critical 
fluctuations themselves and other intersite fluctuations 
are ignored, ^^ SC critical temperatures Tc are given by 




tanh 



^(k) 



2fcBTc 



= 0, (3.24) 



with 



'7r(k) 



_ J cos(fca;a) -I- cos(fcj,a), F = s 
cos(fc2;a) — cos(fcj,a), F = ^7 ' 



(3.25) 



being form factors of the s wave and the d'j wave. Equa- 
tion (|3.24p gives approximate or at least upper limits of 
Tc's in quasi- two dimensions; Tc = +0 K in two dimen- 
sions because of SC critical fluctuations.^^ In quasi-two 
dimensions, when n ~ 1, Tc of the d'y wave is much higher 
than Tc of the s wave since |77d-y(k)p ^ [^^^(k)^ for al- 
most k's except for \kx\ ~ \ky\ on the FS.-^ When I^ < 0, 
triplet p waves are possible: Px and py waves. 

The homogeneous CBW is simply the Fock-type term 
studied in Sec. lIIIBl Since it breaks no symmetry at least 



in the FL, it is not an ordered parameter. When /* is 
weak, the FL is stable against any bond wave. When I^ 
is strong, the FL is unstable against CBW for /^ < and 
SBW for I^ > 0. A flux state, which is a multi-Q bond 
wave with different phases for different Q components, is 
also possible, with Q ordering wave numbers;^ 

When U/W < 1, the perturbation in U is more use- 
ful than that in /s(?w/,q) is. When the nesting of the 
FS is sharp enough, the FL is unstable against an SDW 
state. When U'^Xsi'i'^hfl) is considered as an interaction 
between electrons, it is unstable against an anisotropic 
SC state^^ at least, if no disorder exists. 

The above analysis can never exclude a possibility of an 
exotic state. If it is characterized by an order parameter, 
it is straightforward to study the instability of the FL 
against it. 



IV. DISCUSSION 
A. 1/d expansion theory 

Every term or quantity is classified according to the 
order in l/d in the sitc^^ and wave-numbe r^^i^^ represen- 
tations. The single-site T,„{ie) is of leading order in l/d 
or 0(1). When U/\t\ > Or^ the single-site Xs(iwz) is a 
relevant 0(1) term. Multi-site or intersite terms can also 
be 0(1) only for particular Q's in the wave-number rep- 
resentation. Relevant 0(1) terms are the magnetic Js(Q) 
and Jq^iLUi, Q) for the particular Q's. On the other hand, 
AI]„(ie„,k) = 0(1/ Vd) for any k, Js(q) = 0{l/Vd) and 
jQ{iijJi,q) = 0(l/va) for almost all q's except for the 
particular Q's, and A{iuji,q) = 0{l/d) for any q. 

When the Neel temperature Tn is nonzero for one 
of the particular Q's, magnetization ni(Q) appears at 
T < Tn- Magnetic Weiss mean fields, Js (Q)ni(Q) and 
Jq (0,Q)m(Q), are 0(1). When they are considered in 
the mean-field approximation beyond S'^A or DMFT, the 
mean- field theory is rigorous for d — -l-oo. When T > Tn, 
S^A or DMFT is rigorous for d = -l-oo, except for the 
magnetic susceptibility with the particular Q'S)^ 

The cluster DMFT (CDMFT) is a non-perturbative 
theory for a cluster to include multi-site terms beyond 
S^A or DMFT,Mi^i^i22 In CDMFT, the translational 
symmetry is broken by choosing of a particular cluster 
even in the presence of no order parameter; the symmetry 
is recovered for an infinitely large cluster. The Kondo- 
lattice theory is a perturbative theory to include multi- 
site terms starting from the FL in S'^A or DMFT, which 
is a non-perturbative theory. The conventional pertur- 
bation can treat higher-order terms in l/d; it can also 
treat instability of the FL, as is examined in Sec. IIIICI 
The anomalous perturbation that assumes the existence 
of an order parameter can treat magnetic order, which is 
0(1), and other types of order such as anisotropic super- 
conductivity and bond wave, which are of higher order 
in l/d.'^^ In the Kondo-lattice theory, the translational 
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symmetry is not broken by its framework itself; it is bro- 
ken when an order parameter with nonzero wave-number 
appears. 



B. Magnetism crossover 

The study in this paper is almost restricted to T = OK, 
except for the study in Sec. IIII CI It is straightforward 
to extend the study to T > K. In the Kondo effect, the 
T dependent crossover occurs between a local-moment 
magnet at T > Tk and the FL at T < Tk.'^ In a mag- 
net, there exists a temperature scale Tj^, below which 
magnetic critical fluctuations develop. When Tj^ ^ Tk, 
the magnet is characterized as a typical local-moment 
one. When T^ <C Tk, it is characterized as a typi- 
cal itinerant-electron one. The magnetism crossover is 
simply a crossover between local-moment magnetism at 
T > Tk and itinerant-electron magnetism at T < Tk- 

The local susceptibility Xs(0) gives the Curie- Weiss 
(CW) law for any q or the CW law of a local-moment 
magnet at T > Tk-^^ The RPA polarization function 
Jq(0, q) or P(0, q) gives the CW law of an itinerant elec- 
tron magnet at T < Tk- When there is a sharp nesting 
of FS, P(0, q) gives the CW law only around the nesting 
wave number.— When the chemical potential lies around 
a sharp peak of DOS, P(0, q) gives the CW law only for 
q ~ 0^ The two mechanisms are 0(1) in 1/d. The par- 
ticular q dependences of the two CW mechanisms are 
remnants of those in infinite dimensions. 

When there is no nesting of FS and DOS is almost con- 
stant around the chemical potential, the mode-mode cou- 
pling term — 4A(0, q) can give a local-moment type CW 
law of an itinerant electron liquid at T < Tk.~ When 
there is a sharp nesting of FS or the chemical potential 
is around a sharp peak of DOS, the mode-mode coupling 
term shows an opposite temperature dependence to the 
CW law.s^'™ Since -4A(0,q) = 0{l/d), this mechanism 
does not work in infinite dimensions. 

The magnetism crossover in the periodic Anderson 
model (PAM) is slightly different from that in the Hub- 
bard model. In PAM, conduction electrons and the so 
called d 01 f electrons are strongly hybridized at T <Tk 
to form heavy electrons or quasi-particles, while they are 
independent degrees of freedom at T ^ Tk such that d or 
/ electrons behave as localized spins but conduction elec- 
trons are itinerant. At T <Tk, there is no essential dif- 
ference between the two models. The exchange interac- 
tion Jq [iLUi , q) , which arises from the virtual exchange of 
a pair excitation of quasi-particles, works between quasi- 
particles themselves in either model. At T ^ Tk , there is 
a crucial difference between the two models. In PAM, the 
Ruderman-Kittel-Kasuya-Yosida (RKKY) exchange in- 
teraction Jrkky (iw/ , q) arises from the virtual exchange 
of a pair excitation of conduction electrons, and it works 
between localized spins of d or / electrons in PAM. The 
d or f electron component of the susceptibility of PAM 



is given by Eq. (|3.4p with 

Is{iuJi,q) = Js(q) + Jrkky {ii^h<i) -4A(iw/,q), (4.1) 

at T ^ Tk; the superexchange interaction Js(q) also 
exists in PAM, in general. When T < Tk, Jrkky {ioJi, q) 
turns out to be jQ(iu!i,q).^^ The competition between 
the Kondo effect and the RKKY exchange interaction is 
only relevant at T > Tk or T > Tk. 

C. Evidence for the existence of jQ{iuJi,q) 

The exchange interaction jQ(ia;;,q) or P(?w/,q) is re- 
sponsible for magnetic properties that are observed in 
itinerant-electron magnets. It is similar to the conven- 
tional RPA polarization function except for the pre-factor 
of[4/xl{0)][(j)s/(j)j) . Becauseof this factor, its strength 
is proportional to W* ~ 4fcBTK.— "^ In particular, 

hm P(0, q) = -^ UV4>^) p{0) (X /cbTk, (4.2) 

for the static q = component. Here, Xs(0) — l/(fcB2k) 
is used. Provided that |jQ(0,q)| > |Js(q) - 4A(0,q)|, 
low-energy phenomena are characterized by a single en- 
ergy scale of fcePK, i-e., physical properties obeys the so 
called one-parameter scaling. The one-parameter scaling 
in fceTk is actually observed in the metamagnetic transi- 
tion or crossover in CeRu2Si2r^'2^'^'2& which is evidence 
that jQ(iw;,q) is relevant at least in CeRu2Si2. 

It is obvious that P(0, q) has similar features to those 
of the conventional RPA polarization function. For ex- 
ample, it is antiferromagnetic when the nesting of FS 
of quasi-particles is sharp or the chemical potential lies 
around the center of the quasi-particle band or n ~ 1. It 
is ferromagnetic when the chemical potential lies around 
the top or bottom of the quasi-particle band, i.e., for 
n ~ 2 or n ~ 0. In particular, it is strongly ferromagnetic 
when DOS has a sharp peak at one of the band edges 
where the chemical potential liesr^i^i^^ as DOS's of 
many itinerant-electron ferromagnets such as Fe, Ni and 
so on have. This is consistent with that of Kanamori's 
theory for itinerant-electron ferromagnetism,^^ Since the 
superexchange interaction is antiferromagnetic, however, 
it seems to be difficult for itinerant-electron ferromag- 
netism to occur in the single-band Hubbard model. 

According to Ref. 53, the superexchange interaction 
is ferromagnetic in a multi-band Hubbard model if the 
Hund coupling is strong enough and the band degeneracy 
is large enough. Itinerant-electron ferromagnetism can 
easily occur when both of the superexchange interaction 
and Jq (0,q) are ferromagnetic. 



D. Impossibility of the Mott insulator for finite U 



Since the proof in Sec. IIIBI is made for A — ±0, it 
is another issue what is the ground state in the grand 



canonical ensemble with A = or the canonical ensem- 
ble, where Af or N is & conserved quantity. When A = 0, 
the FS condition may be or may not be satisfied so that 
the ground state of the Anderson model may be or may 
not be a singlet. If it is not a singlet, the ground state 
of the Hubbard model under S'^A or DMFT is an insu- 
lator whose entropy is diverging in the thermodynamic 
limit. In this paper, only such an abnormal insulator is 
called the Mott insulator; if the ground-state entropy of 
an insulator with n ~ 1 or A'^ ~ TV,, is zero or finite in 
the thermodynamic limit, the insulator is called a spin 
liquid. Since it is quite unlikely that the third law of ther- 
modynamics is broken in a relevant Hamiltonian such as 
the Hubbard model with finite U/\t\, we speculate that 
the ground state under S'^A or DMFT is a singlet for 
even A^ or a doublet for odd A^ . Even if the ground state 
can be infinitely degenerate under S'^A or DMFT, the 
degeneracy must be lifted when the RVB mechanism is 
considered beyond it. We also speculate that when U/\t\ 
is finite the ground state within the restricted Hilbert 
subspace must be the FL, an exotic metal, or a spin liq- 
uid for any d in either case of A = and A = ±0. 

According to Lieb and Wu's Bethe-ansatz solution 
for one dimension in the canonical ensemble f^^ when 
N ^ Nf., an M-I transition occurs a,t U — 0. The 
singularity at [/ = is exotic^ The M-I transition 
at [/ = is never due to the opening of the Hubbard 
gap. The insulator for finite U/\t\ is Lieb and Wu's in- 
sulator or spin liquid rather then the Mott insulator. ®° 
When N ^ Nc and U/\t\ > 0, the ground state is the 



Tomonaga-Luttinger liquid 



62.63.64 : 



in which the charge- 



spin separation occurs. When U/\t\ = -|-cxd, in particu- 
lar, the liquid is an abnormal metal whose ground-state 
entropy is diverging in the thermodynamic limit. Be- 
cause of the complete exclusion of double occupancy and 
only the transfer integral t between nearest neighbors 
being nonzero, the charge-spin separation is complete 
such that any eigen-function is a Cartesian product of 
charge and spin parts, i.e., a Slater determinant of spin- 
less fermions and an eigen-function of non-interacting lo- 
calized S = 1/2 spins. ^■^'^■^ The ground-state entropy is 
NkB In 2 and the spin susceptibility obeys the Curie law. 
The charge susceptibility is nonzero and it diverges as 
n ^ 1 - and n ^ -hO at r = K, with n = N/N^. 
When N y^ Nc and U/\t\ — -fcx), spins are localized but 
charges are itinerant. When N =/= N^, the ground state is 
an exotic or abnormal metal. 

In general, electrons are more itinerant in d > 2 di- 
mensions than they arc in one dimension. It is likely 
that the ground state within the restricted Hilbert sub- 
space is a metal at least for 7i ^ 1 or A^ ^^ A^c- In gen- 
eral, the nature of electron correlations is less abnormal 
in d > 2 diinensions than it is in one dimension. It is 
likely that, provided that U/\t\ is finite, the ground state 
within the restricted Hilbert subspace is a singlet in the 
grand canonical ensemble and is a singlet or a doublet in 
the canonical ensemble. 

When the Mott M-I transition is studied in the 



Gutzwiller approximation, the canonical ensemble is con- 
ventionally assumed. The abnormal insulator for N — Nc 
and U > Ubr according to Brinkman and Rice's theory^^ 
must be unstable when the RVB mechanism is consid- 
ered. According to the proof in Sec. Ill Bl any insulator is 
unstable for finite U/\t\ in the grand canonical ensemble 
with A = ±0 even if the RVB mechanism is not con- 
sidered, i.e., under S'^A or DMFT, which is beyond the 
Gutzwiller approximation. 

In numerical S^A or DMFT^i^iH CDMFT^i^SiSe^ 
and Monte Carlo theory j^^^i^^i^^ an M-I transition seems 
to occur when n ~ 1 and U ~ Ubr- The M-I transi- 
tion, at least, for n ^ 1 seems to inconsistent with the 
analysis of this paper, although the phase diagram for 
T > K, which is studied in numerical theories, may be 
different from that for T = K, in general. First of all, 
the effect of A = ±0 or A^ being a non-conserved quantity 
in the grand canonical ensemble, if it is not considered, 
should be explicitly considered also in numerical theories 
because single-particle excitations in a system where A^ 
is not a conserved quantity are different from those in a 
system where A^ is a conserved quantity, as is studied in 
Appendix [C] When an M-I transition occurs in either of 
the numerical theories, a lower-temperature phase seems 
to be the Mott insulator^ It is interesting to examine 
which is actually observed in numerical theories, evidence 
that the third law of thermodynamics holds or evidence 
that it does not. In this context, it is interesting to carry 
out numerical processes beyond S'^A or DMFT in a pa- 
rameter region where the RVB mechanism is expected to 
be effective. If a stabilization effect is observed, it is evi- 
dence that the ground state is a singlet under and beyond 
S'^A or DFFT. If it is not observed, the ground state may 
be infinitely degenerate in the thermodynamic limit, but 
within the restricted Hilbert subspace; the true ground 
state must be an ordered state in the whole Hilbert space. 
If this is the truth, it should be clarified what impedes the 
RVB mechanism or what stabilizes the abnormal ground 
state, where the third law of thermodynamics is broken. 



E. Normal state of cuprate superconductors 

If the RVB state^''-''*' is characterized by an order pa- 
rameter, it is straightforward to study the instability of 
the FL against it. However, no order parameter is pro- 
posed so far. In the mean-field RVB theory;^ the RVB 
state is stabilized by local correlations, which are treated 
with the Gutzwiller projection operator, and the Fock- 
type exchange interaction. No order parameter is intro- 
duced by the Gutzwiller projection operator, which is 
also used to treat the FL, or the Fock-type term, which 
is also nonzero in the FL. On the other hand, it is specu- 
lated that the charge-spin separation occurs in the RVB 
state; bosonic and fermionic elementary excitations are 
called holons and spinous fiS In order to support this spec- 
ulation, a slave-boson RVB theory is proposed in the 
slave-boson t-J model;^'— an electron in the t-J model 
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corresponds to a pair excitation of a slave boson and a 
fermion in the slave-boson t-J model. In this theory, the 
RVB state is characterized by the condensation of slave 
bosons, which are called holons; itinerant fermions, which 
are called spinous^ appear when slave bosons are con- 
densed. However, it is never shown so far what symmetry 
is broken in the t-J model when the condensation of slave 
bosons occurs in the slave-boson t-J model. There is no 
evidence that the RVB state in the t-J model is char- 
acterized by an order parameter, which implies that the 
RVB state, which is an insulator or a metal, is of the same 
symmetry as the FL is. On the basis of the FL theory 
for the Heisenberg model in Appendix [B] which implies 
that the adiabatic continuity^ holds between a spin liq- 
uid in the Heisenberg model and the FL in the Hubbard 
model, we propose that the insulating RVB state is the 
spin liquid and the metallic RVB state is the FL. 

Under S^A or DMFT, p(0) ~ l/W, as is shown in 
Eq. (11321). Beyond S^A or DMFT, p(0) is reduced by 
the RVB mechanism such that p(0) ~ 1/[M^ -f 0-y|cj J|] , 
as is shown in Eq. (|3.18p . On the other hand, DOS's of 
LHB and UHB are not reduced by the RVB mechanism. 
In the so called Hubbard HI approximation^ the band- 
width of UHB and LHB is Whb ^ 0.8W^ for n ~ 1. 
DOS's of LHB and UHB are given by 



p(e, - /i) - P{^a + U-y)~ 1/(2M^hb). 



(4.3) 



When n ~ 1 and n < 1, LHB is just below the quasi- 
particle band. The band centers of LHB and UHB are 
ea- fi^ -Whb/2 and ea + C/ - /i ~ [/ - Whb/^. LHB 
spreads over — Whb ^ e ^ — W*/2, the quasi-particle 
band over |£| < W* /2, and UHB over U-Whb <e<U. 
When UHB is totally above the quasi-particle band, an 
energy region of e ~ —W* is within LHB and that of 
£ ~ W* is within a gap-region between the quasi-particle 
band and UHB, so that p{e ~ -W*) > p{e ~ W*). If 
the reduction of p(0) is small or p(0) ~ 1/M^, 



p{e ~ -W*) ~ p(0) > pie ~ W*). 



(4.4a) 



If the reduction of p(0) is large or p(0) <C l/W, 

p{ec^-W*):$> p{0)c^p{ec^W*). (4.4b) 



An asymmetric p{e) that is consistent with Eq. 
is observed by tunneling spectroscopy;^ Since such an 
asymmetry can arise from the reduction of p(s ~ 0), it is 
evidence that the RVB mechanism is crucial in cuprate 
superconductors. 

In two dimensions, critical fluctuations make Tc down 
to -|-0 K.^'"^ This implies that large deviations from the 
typical FL can occur in anisotropic quasi-two dimensions 
and they are responsible for some of exotic properties 
of cuprate superconductors. If the anisotropy is large 
enough but Tc is still high enough, SC critical fluctuations 
can cause the opening of a pseudo-gap in the SC critical 
regioni21ii2^ It is certain that other effects, such as the 
electron-phonon interaction^'^'^'* and so on, are necessary 
to explain the whole exotic propertiesi^i^ii^i^i^ 



Within the Hubbard model, Tc is low since \J\ is about 
a half of 4|ip/C/ in an actual situation^ Experimentally, 
J is as large as J ~ —0.15 eV.~ In order to explain ob- 
served Tc, the phcnomcnological J ~ —0.15 eV should 
be used in Eq. (|3.7[) or the t-J model with J ~ —0.15 eV 
should be used. The d-p model with relevant parame- 
ters corresponds to the t-J modelM- If Tc should be ex- 
plained microscopically, the d-p model should be used.— 
It is straightforward to extend the study of this paper to 
the t-J model and the d-p model. It is also straightfor- 
ward to extend the study of this paper further such that 
effects of SC, magnetic, and bond wave fluctuations can 
be included, not only in the Hubbard model but also in 
the t-J model and the d-p model. 

According to an early FL theory of high-Tc supercon- 
ductivity in 1987,^^'^^ the normal state is the FL and the 
condensation of d'j-wave Cooper pairs is responsible for 
high-Tc superconductivity. The analysis on possible SC 
states in this paper is simply an extension of the early 
FL theory, or it confirms the early FL theory although 
the extensions are necessary. 



V. CONCLUSION 

The Hubbard model is studied by the Kondo-lattice 
theory. The supreme single-site approximation (S^A), 
which considers all the single-site terms, is reduced to 
determining and solving self-consistently the Anderson 
model, which is an effective Hamiltonian for the Kondo 
effect. It is proved that the ground state under S^A is 
the Fermi liquid except for n = 1 and U/W = -l-oo, 
with n being the electron density per unit cell, U the on- 
site repulsion, and W the band-width. Multi-site terms 
are perturbatively considered beyond S'^A by the Kondo- 
lattice theory. When n ~ 1 and U/W > 1, in particular, 
the Fermi liquid is stabilized under S^A by the Kondo 
effect and is further stabilized beyond S^A by the res- 
onating valence bond (RVB) mechanism. The Fermi liq- 
uid is a relevant normal state to study possible lower- 
temperature phases or the true ground state. In one di- 
mension, the Fermi liquid is unstable against, at least, 
the Tomonaga-Luttinger liquid. In two dimensions and 
higher, the Fermi liquid is unstable against, at least, an 
antiferromagnetic or anisotropic superconducting state. 

It is proposed that the Fermi liquid stabilized by the 
Kondo effect and the RVB mechanism is the normal state 
of cuprate superconductors. In order to explain high 
superconducting critical temperatures, however, the t-J 
model with J ~ —0.15 eV or the d-p model with rele- 
vant parameters, which corresponds to to the t-J model, 
should be used instead of the Hubbard model. 
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APPENDIX A: PRO OF OF THE 
INEQUALITY (f236]l 

Define the following real functions: 

Si{e,e') ^ e + ^I~s' -Re[f:^{s+iO)-iX'^T{s+iO)], (Al) 



S2ie) = -Im[S^(e + iO) - i\^r{e + zO)] , (A2) 



and 



/an I 



^i" (£,£') 



) + Sl{e) 



(A3) 



In this Appendix, it is only assumed that T,a (e + iO) is 
an analytical function in the upper half plane; it may be 
convergent or divergent on the real axis. It follows from 
Eq. ([Ml that 



Y2{e)^l-Sl{e)Y^{e). (A4) 

Since R„ (e + iO) defined by Eq. (|2.25p is given by 

R„{e + iQ)^Yi{e)-iS2{e)Yo{e), (A5) 

the mapping condition (j2.22bp is given by 

Y^{e)-Yo{e)Y2{e) 



A(e) = ReA^r(e+iO) - 52(e) 



Y^{e) + Sl{e)Yi{e) 



(A6) 



Here, Eq. (|A4p is used. It is trivial that 5*2 (e) > and 
Y^{£) + Sl{e)YQ{e) > 0. Since an inequality of 



.,. [x + S^{e,e')Y 



h'^^'^U^ 



) + Sl{e) 



> 0, (A7) 



or Yo{e)x'^ + 2Yi{e)x + ^2(2) > is satisfied for any real 
X, it follows that 

Y^{e)-Yo{e)Y2{e)<0. (A8) 

The inequality (P?^ . A(£) > ReA2r(e + iO), holds. 



APPENDIX B: FERMI-LIQUID THEORY FOR A 
SPIN LIQUID IN THE HEISENBERG MODEL 

1. Localized spin in s-d model 

The FL theory for the s-d model, which is a prototype 
of the FL theory for the Heisenberg model, is first studied 
by a different approach from Nozieres'.— In the s-d limit, 
which is defined by f//|y|2 -> +00 with 2{ed- fl)-\-U = 
and Js-d — — 2|yp/[/ kept constant, the Anderson model 
([2?T4l) . but with constant V{k) = V and 



^^<5[£ + /i-i?,(k)]=i5,(0). 



Na 



(Bl) 



is reduced to the s-d model: 
Hs-d = -Js-d^(S • cr'''' )cJ^cot 



X!^c(k)4^Ck^, 



kCT 



(B2) 

with S a localized spin with S = 1/2 at the 0th site and 
cqt — (1/v^) Sk'^kr- A constant term of Ea is ignored 
in Eq. (|B2|. 

It is straightforward to extend the FL theory for the 
Anderson model''' ^ to the s-d limit. The Green functions 
for d electrons and conduction electrons in the Anderson 
model are given, respectively, by 



G^(e + zO) = -3- 

07 



1 



iA* 



0{s^ 



(B3) 



and 



Gkk'a{e-\-iO) = gka{£+iO) + 9kai£+iO)gk'ai£+iO) 



A* 



■^i^c(O) 



1 



e + iA 



7 + Oie') 



for |e| < fceTk, with 



A* = ^n\V\^D,{0) = -kBTK, 



(B4) 



(B5) 



and 5k<7(£ + «0) == l/[e + A - -Bc(k) + iQ]. In Eq. (|B5|) . 
Eq. (|2.29p is used. Since Tk is nonzero and finite in the 
s-d limit, 

^-y oc \V\'^ oc U, (B6) 

is satisfied there. In the s-d limit or in the limit of 0^ — > 
+CXD, the Green function of d electrons is vanishing, as is 
shown in Eq. (jB3[) , but the fermionic spectrum of single- 
particle excitations for d electrons is still well defined, as 
is shown in Eq. (jB4p . 

In the s-d model, d electrons are exactly localized and 
they carry a localized spin. It is trivial that the single- 
particle or fermionic Green function of the localized spin 
can never be defined. However, the fermionic spectrum, 
which describes scatterings of conduction electrons by the 
localized spin, is defined by Eq. (|B4p . Physical properties 
of the localized spin in the s-d model can be described 
by the fermionic spectrum according to the FL relation. 
For example, the specific heat is given by C = ^T -)-••• 
at T < fee Ik, with 7 given by Eq. (H^D]). This is simply 
the FL theory for the s-d model by Nozieres.^^ 

The number of d electrons that carry a localized spin 
is a conserved quantity in the s-d model but the number 
of d electrons is not in the Anderson model, i.e., local 
gauge symmetry exists in the s-d model but it does not 
in the Anderson model. It should be noted that, how- 
ever, the adiabatic continuity^ holds between the FL in 
the s-d model and the FL in the Anderson modeli^ It 
should also be noted that the Kondo peak, to which the 
Gutzwiller band in the Hubbard model under S "^A or 
DMFT corresponds, appears in the fermionic spectrum 
even if no d electron can be added or removed in the s-d 
model but when an electron can be added or removed in 
the conduction band, as is shown by Eq. (|B4p . 
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2. Spin liquid in the Heisenberg model 

In the Heisenberg limit, which is defined by U/\t\ -^ 
+00 with 2(ea - ^) + [/ = and J = -At^ /U kept con- 
stant, the Hubbard model with the electron reservoir, 
which is defined by Eq. (12. 1|) . is reduced to the Heisen- 
berg model with a thermal reservoir: 



(y> 



-^'^'■EE(s^"^''>l^-' 



(B7) 



ieTZ tt' 



with Si a localized spin with S = 1/2 at the ith site 
and Jr = —2\v\'^/U. A constant term of Nc{ea — /i) is 
ignored in Eq. (|B7p . First, two dimensions and higher are 
assumed. When the RVB mechanism is only considered 
beyond S^A or DMFT, the ground state is the FL even 
in the Heisenberg limit provided that the perturbation 
X^Jr from the thermal reservoir is weak enough. Then, 
it is assumed that 



and 



lim |wp/C/ < +00. 

f7/|t|-»+oo 



lim \v\'^ /(py < -fcxD, 

U/\t\^+oo 



(B8a) 



(B8b) 



with (jj-y one of the expansion coefficients for the self- 
energy renornralized by the Fock-type term; the model 
(|B7p and quasi-particles are well defined when Eq. (IB8|) 
is satisfied. It is also assumed that the FL is totally stabi- 
lized by the RVB mechanism, i.e., it is assumed according 
to Eq. (|3.19p and W oc \t\ ex ^/U that (f)j satisfies 



(B9) 



lim VU/<j)y = 0. 

(7/|iH+oo 



The Green functions for electrons in the Hubbard 
model and the reservoir averaged over the ensemble are 
given, respectively, by 



and 



1 



1 



^(k)+iX2T*{e+iO) 



0(£2 



(BIO) 



Gbais+iO, k) = gb„{e + iO, k) + nhA^^g^^(e + iO, k) 

1 



s-S.O<.) + iX^r*{e+iO) 



-0{e' 



, (BU) 



for |e| < |cjJ|, with_£b^(i£„,k) = l/[ie„ + ^l - £'b(k)], 
f (k) given by Eq. ([OI]) . and 

r*(£ + ^0) - jnt^-^ V5fc,(e -I- iO, k). (B12) 



In the Heisenberg limit, the Green function for electrons 
is vanishing. Provided that Eq. (|B8b|) is satisfied, how- 
ever, quasi-particle excitations are well defined. 

The conductivity in the Heisenberg limit is given by 



e 

Nr. 



d_ 
dkx 



<J.M oc _^_iJ(k)^[i?(k)-fAI].(k)] 



d_ 
dkx 



\t\ 



|t;|2A2 



X [lmG^{+iO,-k)Y 
ne^\Jl 



(B13) 



where the vertex correction consistent with the Fock-type 
self-energy is included. The asymptotic behavior of Jr 
or I up should be properly assumed to satisfy Eq. (jB8[) . 



Assume that (/)-y (x U^, where k > 1/2 is required by 
Eq. jHll. If 1/2 < K < 1, \v\^ ex t/'* is assumed; r*(£+iO) 
is nonzero and finite and Jr is vanishing. If k = 1, jwp oc 
U is assumed; both of r*(e -I- iO) and Jr are nonzero 
and finite. If k > 1, jup ex [/ is assumed; r*(£ -I- iO) is 
vanishing and Jr is nonzero and finite. In either case, the 
conductivity is vanishing because jil/diipA^) -^ in the 
Heisenberg limit followed by the limit of A -^ 0. 

Assume that A^ is small but nonzero. When 1/2 < 
K < 1, X^r*{e + iO) is nonzero for vanishing X^Jr, which 
means that there is a singularity at A^J,, = 0. When 
|u|2 (X U is assumed, for example, X^r*{e + iO) is diverg- 
ing for nonzero and finite X^Jr, even if X^Jr is infinites- 
imally small. The ground state is completely disordered 
due to an infinitesimally small X^Jr- This is unreason- 
able. When K > 1, }?T*{£ + iO) is vanishing for nonzero 
and finite X^Jr, even if }?Jr is large. A finitely large per- 
turbation has no effect on the ground state. This is also 
unreasonable. We speculate that the truth is k = 1 or 
(j)y oc [/, as it is in the s-d limit for the Anderson model. 

When the RVB mechanism is only considered beyond 
S'^A or DMFT within the restricted Hilbert subspace, 
the ground state in the Heisenberg limit is a spin liq- 
uid, where the Green function of electrons and the con- 
ductivity are vanishing but the fermionic spectrum ^(k) 
of almost localized electrons is defined. It is probable 
that, under the corresponding approximation within the 
restricted Hilbert subspace, the spin liquid is also the 
ground state of the Heisenberg model, where the single- 
particle Green function of localized spins can never be 
defined and the conductivity exactly vanishes. Although 
^(k) related to localized spins is also defined, it does not 
exist in the Heisenberg model. Physical properties of the 
spin liquid can be described by ^(k) according to the FL 
relation. ^^J^ For example, the specific heat per unit cell 
is given by C = ^T + ■ ■ ■ at T <C |c,/ J|, with 7 given by 
Eq. (|2.30p . This is simply the FL theory for the Heisen- 
berg model. 

The FS exists in the Hubbard model but it does not in 
the Heisenberg model. In general, no change of symmetry 
occurs between a metal, in which the FS exists, and an 
insulator, in which no FS exists, if no order parameter 
is involved; the breaking of local gauge symmetry may 
occur/SS. The analysis in this Appendix implies that the 
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adiabatic continuity^^ holds between the FL in the Hub- 
bard model and the spin liquid in the Heisenberg model, 
i.e., the spin liquid is simply the FL. We propose that the 
insulating RVB stated is the spin liquid and, therefore, 
the metallic RVB statei^ is simply the FL. 

In one dimension, terms proportional to elnjej ap- 
pear in the multi-site self-energy at T = K. Even in 
this case, it is possible to describe the Green functions 
by Eqs. (|BT0)) and (|BTT|) . where each of l/[e - C(k) + 
i\'^r*{e + iO)] is replaced by l/[e + cJe\n\e\ - ^(k) + 
iX^r*{e+iO)], with c a numerical constant. It is straight- 
forward to extend the above analysis to one dimen- 
sion. The extended analysis implies that the adiabatic 
continuity^ holds between the Tomonaga-Luttinger liq- 
uid in the Hubbard model and Bonner and Fishers's spin 
liquid^"" in the Heisenberg model. 

The Tomonaga-Luttinger liquid and Lieb and Wu's 
spin liquid seem to be of the same symmetry as each 
other. It is interesting to examine whether the adiabatic 
continuity2^ holds between them. If the adiabatic conti- 
nuity holds between the Tomonaga-Luttinger liquid and 
Bonner and Fishers's spin liquid, as is discussed above, 
it is probable that the adiabatic continuity also holds 
between the Tomonaga-Luttinger liquid and Lieb and 
Wu's spin liquid. When A = ±0, A^ is a non-conserved 
quantity. When U/\t\ is finite and the chemical poten- 
tial is continuos at n = 1 as a function of n, a state for 
71 = 1 is, in a certain sense, an average over states with 
N = Nc + AN, with AiV = 0, ±1, ±2, • • • , in the thermo- 
dynamic limit. One may argue that if the ground states 
for n ^ 1 ± are the same as each other the ground 
state for n = 1 must be the same as those for n ^ 1 ± 0. 
If this argument is relevant and the adiabatic continuity 
holds between the Tomonaga-Luttinger liquid and Lieb 
and Wu's spin liquid, it is possible that the ground state 
for A = ±0, n = 1, and finite U/\t\ in one dimension is 
simply the Tomonaga-Luttinger liquid. 



with N the number of total electrons and a a quantum 
number for an eigen-state of Ti.. The ground state for /i 
that corresponds to n = 1 is denoted by |iVg;g). The 
thermodynamic limit of Ng — > -|-oo is assumed; /i's are 
the same as each other among N — Ng, NgzLl, ■ ■ ■ . 

When \t\/U -^ +0, it follows that (^ - ea)/|i| ^ +oo 
and {ca + U — fi)/\t\ — + -|-oo. All the unit cells in the Hub- 
bard model are singly occupied in either of the ground 
states for N ~ Ng, Ng ±1, • • • ; extra electrons go to 
the reservoir or deficient ones come from it. They are 
degenerate with each other such that Ej^ = Ejq ^^ 
— ■ ■ ■ , which is simply denoted by Eg. Eigen-states are 
classified according to the numbers of doubly occupied 
and empty unit cells in the Hubbard model, which are 
denoted by d^v-a ^'^'^ ^N-ai respectively: 

En-c =Eg+ dj^,^ {ea + U - ii)+ Cf, (^ - ea) , (C2) 



for N = Ng, Ng±l, ■ ■ ■ . They are degenerate with each 
other when their dff.^ and e^.^, are the same as each 
other, even if their iV's are different from each other. 

When < \t\/U < 1, it follows that iJ.-€a ^ U/2 > |i| 
and €a + U — ^ ~ U/2 :§> \t\. Eigenstate can be still 
classified according to dj^.^ and ej^Q. The degeneracy 
(|C2p in the ground and low-lying states with d^va — 
and e^.^ = is lifted by second-order perturbation in t: 



E 



N:i 



Eg±Oi\j\), 



(C3) 



for N — Ng, Ng ± 1, ■ ■ ■ , with J the superexchange in- 
teraction. It should be noted that the ground states for 
N — Ng, Ngztl, ■ ■ ■ are still degenerate with each other. 
The 0(| J|) term in Eq. (jCSp implies that each of the 
ground states is a singlet stabilized by the RVB mecha- 
nism. Since a doubly-occupied or empty site is itinerant, 
the degeneracy (jC2[) in excited states with dpf.^ > 1 or 
^N-a > 1 is lifted by first-order perturbation in t: 



APPENDIX C: SINGLE-PARTICLE OR 
FERMIONIC EXCITATIONS 

1. Grand canonical ensemble 

One of the purposes of this Appendix is to study a 
general feature of single-particle or fermionic excitations 
in the grand canonical and canonical ensembles under 
an assumption that d > l,n~ lorA^~A^c, and 
1 <C U/\t\ < -|-oo. It is shown that fermionic excita- 
tions, which themselves are never observables, are differ- 
ent between the two ensembles. Then, it is shown that 
the well-known physical picture for the Mott transition, 
which is one in the canonical ensemble, is never relevant 
to explain the Mott transition. 

First, consider the Hubbard model in the grand canon- 
ical ensemble or H defined by Eq. (|2.ip . The eigen- 
equation is given by 

(CI) 



n\N;a) =Ej^.^\N;a), 



En;c. ^Eg + d^.^ (ea + U-fi) + e^.„ (^-e,) ± O(W^), 

(C4) 



iOY N = Ng, Ng±l, 



When T = K, DOS defined by Eq. (P?^ is given by 



Pie) = 



ka 



S(e-E, 



E, 



(C5) 



with Nc the number of unit cells of the Hubbard model 



and 



^k« = |(^9 + l;"|<|^.;5)| , (C6a) 



A^a^\{Ng-l;a\a^,\Ng;g)\\ (C6b) 



In the grand canonical ensemble, the fermionic spectrum 
p{e) is related to response functions or observables, which 
are bosonic, according to the FL relation . ^^i^^ 
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Low-energy excitations are possible from the ground or 
initial state \Ng;g') to low-lying final states INg ± l;a) 



with d 



N„±l;a 



and 



-N„±l;a 



0. A narrow band, 



whose width is 0(|J|) according to Eq. (|C3|) . appears 
around e ~ or around the chemical potential. The 
nature of the band depends on that of the ground state. 
When the ground state is the FL, the band is simply the 
Gutzwiller band renormalized by the RVB mechanism. 
A band or structure also appears, or at least |i|/o(0) > 0, 
in the Tomonaga-Luttinger liquid. If the ground state is 
a spin liquid, |t|p(0) — but a structure must appear in 
p{e) at \e\ < 0(| J|). It is also possible that a complete 
gap opens if an order parameter appears. 

Even if the ground state is a spin liquid, the charge- 
spin separation occurs, or an order parameter appears, 
excitations as large as U are possible from the ground 

1. Two broad 

bands, whose widths are 0{W) according to Eqs. (|C3p 
and (|C4[) . appear around e ~ ±C//2 below and above 
the chemical potential. They are simply UHB and LHB. 
When U >W, the Hubbard gap or pseudo-gap opens for 
any d > 1 between UHB and LHB. 



Canonical ensemble 



state to final states with dj^ _^i.^ — 1 and e^ _^i.^ — 
and those with dj^ _i.a — and e^y _i.^ 



for the ground and low-lying states. 

The 0(|J|) term in Eq. (|C9p implies that the ground 
state for iV = A^c is a singlet for even Nc or a doublet 
for odd A^c stabilized by the RVB mechanism. It also 
implies that spin fluctuations with energy scale 0(| J|) are 
developed. When no charge-spin separation occurs, spin 
fluctuations inevitably couple with charge fluctuations. 
The existence of coupled charge-spin fluctuations implies 
that the ground state is a metal. When the charge-spin 
separation occurs, the ground state may be a spin liquid 
or a metal. The ±0{W) term in Eq. (IC11|) . which is 
due to the itineracy of a doubly-occupied or empty site, 
implies that either of the ground states for A^ = A^c ± 1 is 
a metal. A phase diagram speculated from this argument 
is consistent with that for one dimension according to the 
Bethe-ansatz solution, where the ground state for finite 
[/ is a spin liquid for N — Nc but is a metal for N ^ Nc- 
The argument implies that the ground state for finite U 
is a metal for any A^ in two dimensions and higher, if 
no charge-spin separation occurs or no order parameter 
appears. 

Consider the conductivity a^xi^) for A^ = A^c- It is 
mainly determined by two-particle excitations. When a 
pair of single-particle excitations, a particle and a hole, 
are bound, i.e., almost all unit cells are singly occupied 
in a pair-excited state, it follows that 



Next, consider a Hubbard-like model in the canonical 
ensemble: 



Tiu' =n + U' {Af - Ncf 



(C7) 



with H the Hubbard model (|2.2p and JV the number op- 
erator (|2.5[) . Since A/" is a conserved quantity, none of 
observables depend on U' since U' {Af — Ac) is simply 
a constant of U' {N — A"c) . The eigen-equation is given 
by 



nu'\N;a) ^EN;c.\N;a). 



(C8) 



Since almost all the unit cells are singly occupied in the 
ground and low- lying states for N = Nc, it follows that 



EN^,^ = Ncea±0{\J\), 



(C9) 



by second-order perturbation in t. When a pair of empty 
and doubly occupied site appear in an excited state. 



E 



Na 



Ncea + U±0{W), 



(CIO) 



by first-order perturbation in t. When an electron or a 
hole is added, it remains within the Hubbard model, which 
is sharp contrast with that it can escape from the Hubbard 
model in the grand canonical ensemble. By first-order 
perturbation in t. 



SjV. + l:a ^U'+ {Nc + l)ea + U± 0{W), (Clla) 



and 



E 



Nc~l:a 



U' + {Nc-l)ea±0{W), 



(Cllb) 



E 



Na 



E 



N,;g 



0{\J\), 



(C12) 



with the ground state denoted hy a — g. When a particle 
and a hole are not bound in an excited state, 



E 



N, 



E 



Na 



U±OiW). 



(C13) 



None of pair excitations depend on U' . Two struc- 
tures can appear in Keaxxi^)' a low-energy one around 
uj = 0{\J\) and a high-energy one around uj = U . If the 
ground state is a spin liquid, a small gap opens within the 
low-energy structure. The small gap is 0(| J|) or smaller 
than 0(| J|). If the ground state is a metal, no gap opens. 
If the ground state is the FL, in particular, the low-energy 
structure is simply the Drude term due to quasi-particles, 
whose band-width is 0(| J|). The high-energy structure 
corresponds to a pair excitation between the Hubbard 
gap. The analysis on the conductivity, which is an ob- 
servable, in the canonical ensemble is consistent with that 
on single-particle excitations in the grand canonical en- 
semble, which are related to observables according to the 
FL relation,iLi^ 

A fermionic spectrum is defined by 

ka 



+ Kj{£ + EN^-l;C~EN^;g) , (C14) 



with 



■^ka 



\{Nc + \\a\al^\Nc]g)\ 



(C15a) 
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and 



^ka = |(^c-l;a|ak^|iVc;.9)| 



(C15b) 



A gap as large as U + 2U' — 0{W) opens in Xp{e). When 
U' — 0, TCu' is simply the Hubbard model Ti.. One may 
argued that the Hubbard gap as large as [/ — 0{W) or 
U — W opens in Xp{s) so that the ground state for N = 
Nc is the Mott insulator for U > W. This is simply the 
well-known physical argument or picture for the Mott 
transition. However, the picture is never relevant. First 
of all, Xp{e) depends on U', which means that Xp{e) is 
not related to observables or, at least, it is not directly 
related to observables. When the ground state is a spin 
liquid, a gap in Xp(s) is different from the small gap in 



the conductivity Reaxx{i^)- When the ground state is a 
metal, Xp{e) cannot describe observables of the metal, 
whose energy scale is 0(| J|). 

All the analyses and arguments in this paper show or 
imply that the Mott insulator is impossible when U/\t\ is 
finite, J is nonzero, or the RVB mechanism is effective. 
An exception is the well-known physical picture. The 
FL theories for the s-d model and the Heisenberg model 
in Appendix |B] imply that a relevant fcrmionic spectrum, 
which is related to observables, can be defined even in the 
canonical ensemble when a thermal reservoir is explicitly 
considered, i.e., when an electron is added or removed 
in the thermal reservoir even if no electron is added or 
removed in the Hubbard model. 
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